Popular summary of: 

Ozone Loss from Quasi-Conservative Coordinate Mapping during the 
1999 — 2000 SOLVE Campaign 

During the winter of 1999-2000, the Sage III Ozone Loss and Validation 
Experiment (SOLVE) field experiment took place in Kiruna, Sweden. The 
purpose of SOLVE was to examine ozone depletion mechanisms in the Arctic 
stratosphere (from about 10 to 50 km altitude) during the winter and 
early spring, when a band of strong winds (the "polar vortex”) circle 
the pole. Measurements of stratospheric ozone were made by several 
different kinds of instruments in different meteorological situations. 

We analyzed these data using the "quasi-conservative coordinate mapping- 
technique, in which the measurements are analyzed m terms of 
meteorological properties ("potential temperature" and "potential 
vorticity") which tend not to change very much over a few days. This 
technique reduces or removes the changes that are associated with the 
polar vortex moving around. Over longer time periods, potential 
temperature and potential vorticity change as air cools and descends 
within the polar vortex. We account for these changes by calculating 
the trajectories of air parcels, and this enables us to extend the 
analysis over a ten-week period from January 10 to March 17, 2000. Using 

data from the NASA ER-2 aircraft, from the DIAL and AROTEL laser 
sounders on the NASA DC-8 aircraft, and balloon-borne ozonesondes, our 
analysis reveals changes in ozone which, because we have removed the 
effects of polar vortex motion and the descending air, indicate chemical 
destruction of ozone in early 2000. We find a peak decline rate of 
approximately 0.03 ppmv/day near 470 Kelvin of potential temperature 

(near 20 km) in mid-January which sinks in altitude to around 440 Kelvin 

(near 18 km) in mid-March. 
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Abstract. 

Ozone observations from the AROTEL and DIAL lidars on board the NASA DC-8 
aircraft, the NOAA in situ instrument on board the NASA ER-2 aircraft, and THESEO 
2000 ozonesondes are analyzed using a quasi-conservative coordinate mapping technique. 
Measurements from the late-winter/early-spring SOLVE period (January through March 
2000) are incorporated into a time-varying composite field in a potential vorticity-potential 
temperature coordinate space; ozone loss rates are calculated with and without diabatic effects. 
The loss rate from mid-January to mid-March near the 450 K isentropic surface in the polar 
vortex is found to be approximately 0.03 ppmv/day. 
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1. Introduction 

Since the discovery of the Antarctic ozone hole by Farman et al. [1985], stratospheric 
researchers have monitored the ozone layer, watching for signs of similar chemical processes 
over the Arctic regions. To detect ozone loss, researchers have collected and analyzed 
measurements from satellite instruments such as the Total Ozone Mapping Spectrometer 
(TOMS) and Microwave Limb Sounder (MLS). {Newman et al. [1997] and Manney et al. 
[1997] are examples.) Just as important, however, have been data from field experiments 
such as the Airborne Arctic Stratospheric Expeditions (AASE I and II), as outlined in 
Tu rco et al. [1991] and Rodriguez [1993]. Measurements from such experiments typically 
involve a large number of species from relevant families of atmospheric constituents, and 
these data — frequently of high spatial and temporal resolution — enable a fairly detailed 
understanding of the chemical processes associated with stratospheric ozone loss. (See. for 
example, Salawitch et al. [1993].) 

During the winter of 1999-2000, the Sage III Ozone Loss and Validation Experiment 
(SOLVE) was carried out simultaneously with the Third European Stratospheric Experiment 
on Ozone (THESEO) campaign. Meteorological and trace gas constituent data were collected 
by many different in situ and remote sensing instruments. These included measurements 
of stratospheric ozone by airborne lidars, an airborne UV-absorption photometei. and 
ozonesondes. 

Given SOLVE’s and THESEO's variety of ozone measurements taken by different 
instruments at different times and scattered locations, it is desirable to put these data into 



a consistent wide-scale meteorological framework. As outlined in Schoeberl and Lait 
[1991], constituent mapping using quasi-conserved coordinates can be a useful technique 
for this purpose. A composite field of trace gas constituent measurements is constructed in 
a vortex-relative coordinate system using potential vorticity (PV) and potential temperature 
(6) as abscissa and ordinate, respectively. Then, given gridded meteorological fields of PV 
and 9, one can map the composite into real space, creating a three-dimensional field of the 
reconstructed constituent. Schoeberl et al. [1989] and La it et al. [1991] demonstrate examples 
of this “P V-0” technique. Related techniques are employed in Manney et al. [1994] and Kyro 
etal. [2000], 

In this paper, we refine the PV-fl technique and use it to investigate ozone depletion 
observed during the late winter/early spring of 2000. Section 2 describes the data used in this 
analysis, and Section 3 describes the enhanced technique. Results are presented in Section 4, 
and the conclusions follow. 

2. Data 

2.1. Meteorological Data 

PV-0 analysis requires that each measurement be associated with a value of PV and 0. 
and PV at least must be obtained from gridded meteorological analyses. Three sources ot such 
analyses were used for this work. 

The NASA Goddard Space Flight Center’s Data Assimilation Office (DAO) product tor 
the SOLVE period is obtained from their GEOS-3 assimilation system for EOS-Terra suppoit. 
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GEOS-3 is the successor to the GEOS-1 system documented in Pfaendtner et al. [1995]. 
These data grids extend from 1000 hPa to 0.2 hPa, have a horizontal resolution of 1° longitude 
by 1° latitude, and are produced four times daily. 

The United Kingdom Meteorological Office (UKMO) product generated for the Upper 
Atmosphere Research Satellite (UARS) project is another assimilation effort, described in 
Swinbcink and O’ Neil [1994], The version used here extend from 1000 hPa to 0.4 hPa, have a 
horizontal resolution of 3.75° longitude by 2.5° latitude, and are produced once per day. 

The third source of meteorological data is from the long-term assimilation performed by 
the National Centers for Environmental Prediction (NCEP) and the National Center for Atmo- 
spheric Research (NCAR) The NCEP/ 

NCAR reanalysis system’s procedures are applied consistently to over 40 years ot raw 
data, resulting in a dataset which is useful for long-term studies. (For further information, see 
Kalney et al. [1996]). These data extend from 1 000 hPa to 10 hPa, have a horizontal resolution 
of 2.5° longitude by 2.5’ latitude, and are produced tour times daily. 

To account for diabatic effects, heating rates were calculated from the UKMO analyses 
using the model described in Rosenfield et al. [1994], 

2.2. Ozone Data 

SOLVE consisted of three deployments: December 2-14, 1999; January 10-February 3. 
2000; and February 26-March 1 6, 2000. During these periods, the NASA DC-8 and ER-2 
aircraft flew numerous sorties out ot Kiruna, Sweden. 

On board the ER-2, a dual beam UV-absorption photometer measured ozone volume 
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mixing ratios at flight altitudes (above 18 km). Proffitt and McLaughlin [1983] describe the 
instrument. These data are well-calibrated and of good accuracy (around 3%). And at one 
measurement per second on an aircraft moving at approximately 200 m/s, they have high 
spatial and temporal resolution. However, their coverage is limited to the actual position of the 
aircraft. 

In contrast, two lidar instruments on board the DC-8 measured profiles of ozone 
above the aircraft well into the stratosphere. The UV Differential Absorption Lidar (DIAL) 
measurements have a vertical resolution of approximately 750 m and are taken at least five 
minutes apart (giving a horizontal resolution of up to 70 km). DIAL profiles reach to 25 or 30 
km. DIAL measures ozone number densities, and the DAO analyzed temperatures are used 
to convert to volume mixing ratios. See Browell et al. [1998] for further description of the 
instrument. 

The Airborne Raman Ozone and Temperature Lidar (AROTEL) uses Rayleigh and 
Raman scattering to obtain profiles with a vertical resolution of 0.5 to 1.5 km and which 
are taken at least 2 minutes apart (giving a horizontal resolution of up to 24 km). AROTEL 
profiles reach up to 35 to 40 km. Like DIAL, AROTEL measures ozone number densities, but 
AROTEL uses its own temperature measurements to convert to volume mixing ratios. McGee 
et al. [in this issue] describe this instrument in more detail. 

These three aircraft-based instruments provide stratospheric ozone measurements of high 
resolution along their flight tracks which cut cross large areas ot the Arctic between Greenland 
and Novaya Zemlya. They arc limited to flight days within the three SOLVE deployment 


periods, however. 
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Ozonesondes launched by the THESEO campaign, the Canadian Atmospheric 
Environment Service, the World Meteorological Organization network, and Russia measured 
profiles at a number of sites in the Arctic. Using an electrochemical concentration cell, these 
instruments take measurements with a vertical resolution of 8 to 50 m and typically reach up 
to around 30 km. Although the launch sites are geographically sparse, these profiles exist for 
nearly every day of the period examined here. 

3. Analysis 

PV was obtained from the meteorological analysis closest in time to each measurement by 
interpolating the gridded field bihnearly in the horizontal and using cublic spline interpolation 
in log-pressure for the vertical. The lidar instruments profiles are a function ot geometric 
height above the aircraft', the analyzed geopotential heights on pressure sui faces weie 
converted to geometric heights and used to obtain the pressure of each measurement along 
the flight track . Modified PV (Lait [1994]) was used to avoid the strong vertical scaling 
exhibited by regular Enel's PV. For consistency, the 0 values were similarly interpolated from 
the meteorological analyses, even tor instruments which measured temperature. 

A time-varying ozone composite was constructed on a grid in PV-0 space. The edges of 
the grid were chosen to be well-removed trom the values typical of the area of intemst (the 
Arctic lower stratosphere): 0 to 50 PV Units (One PVU being one 10 f, K m 2 /kg - s) and 200 
to 1000 K in 0. The grid coordinates were nondimensionalized to a rectangle with an aspect 
ratio on the order of unity. 

A reasonably long-lived trace gas is expected to be well-mixed along contours of PV on a 
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surface of constant 9 (Leovy et al. [1985]), and in the absence of diabatic affects and chemical 
changes a simple composite made by averaging mixing ratios near a gridpoint should yield an 
accurate picture of a time-invariant trace gas distribution in PV-0 space. 

But over the course of the ten weeks examined here one cannot neglect diabatic effects or 
chemical ozone loss. Both effects appear as a change in ozone mixing ratio at a given point in 
PV-0 space. One must take into account this time-varying nature of the field when creating the 
PV-0 composite. 

Over periods of a few weeks (shorter than seasonal time scales), the ozone mixing ratio \ 
at fixed PV and 9 values can be approximated as a linear change in time t: 

X(P V,0) = a(PY,0) + 6(PV. 9)t (1) 

To make the ozone P V-0 composite, constants «(PY,0) and 6(PY,0) were computed by 
constructing a time series at each gridpoint of all the nearby ozone measurements. A linear 
least-squares fit in time was applied to the measurements, which were weighted by 

»i» i = (l + l/v^)exp(-d?) (2) 


where a t is the estimated or quoted uncertainty in the ith measurement and d t is a scaled 
nondimensional distance from the /th measurement to the gridpoint in P V-9 coordinates. A 
typical time series and fit for the ozonesonde data is shown in Figure I. 


Figure 1 


The results of the time fits are two PV-0 grids a and b , containing slopes and intercepts, 


respectively. An ozone field in PV-0 space can be reconstructed from these for any day. Given 
PV and 9 fields at points in real space, this reconstructed ozone grid can be mapped onto those 
points to obtain nominal values oi ozone there. Figure 2, tor example, shows ozonesonde data 


Figure 2 
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which has been mapped onto vertical profiles above the NASA DC-8 tor comparison with the 
AROTEL data. The figure shows the comparison for the flight of March 9, when the DC-8 
flew back and forth across the vortex edge. The remapped sonde data generally compare well 
with the AROTEL measurements. Two properties of the reconstructed field stand out. First, its 
spatial resolution is much less than that of the original data. This is not surprising, considering 
the averaging process and the limited resolution of the meteorological analyses. Nevertheless, 
the reconstructed field does show the edge structure clearly. Second, the reconstructed 
ozonesonde data cover a much larger part of the flight than the actual AROTEL measuiements. 
This is true despite the fact that the raw sonde profiles are much more sparsely distributed 
over the Arctic region than the AROTEL profiles. The reconstruction process assumes that the 
average of ozone values measured under certain meteorological conditions will be typical of 
all ozone values in identical conditions, and maps them accordingly. 

In a similar way, one can use gridded three-dimensional fields of PV and 0 to map the 
ozone composite onto the fields' grid points to obtain an ozone field covering a large area. 
After using climatological values to fill in regions above and below the reconstructed profiles, 
such a field can be integrated vertically, yielding a total ozone field that can be compared 
with TOMS measurements. A comparison for a typical day is shown in Figure 3. PV-0 
reconstruction reproduces the morphology of the total ozone field fairly well, and the values 
themselves are reasonable. The lower total ozone values from the reconstruction are higher 
than the TOMS measured values. Inspection of ozonesonde and reconstructed ozone profiles 
indicates that these differences are chiefly caused by uncertainties in the reconstructed values 
in the lowermost stratosphere/upper troposphere, where small variations in ozone mixing 


Figure 3 



ration can result in large variations in total ozone. 


Reconstructed fields can be used to determine ozone change rates by differencing 
the fields at the beginning and end of a time period, or by mapping the slopes from the 
least-squares fits into real space. However, the change caused by diabatic descent of air in the 
polar vortex and the change caused by chemical destruction of ozone will be indistinguishable. 
A way is needed to account for the diabatic descent to enable the rate of chemical destruction 
to be determined. Two methods for accounting for diabatic effects were used in this analysis. 

To determine how much descent took place within the vortex over this period, air parcel 
trajectories were traced using an isentropic trajectory model with diabatic corrections (See 
Schoeberl et al. [1998], but here the net diabatic heating rate is not balanced). Parcels were 
initialized at approximately equally-spaced gridpoints on 6 surfaces spaced 10 K apart from 
400 K to 600 K, and the trajectories were run both forward and backward between January 
10 and March 18 using the UKMO meteorological data and heating rates from the Rosenfield 

model. 

To obtain descent rates within the vortex core, only those parcel trajectories which 
started out and ended up in the core of the vortex were retained. The vortex core was found 
by determining the PV value that marks the highest 10% all PV values on the (equidistant) 
gridpoints which lie poleward of 50° N on each theta surface. These PV values were used to 
delimit the core of the vortex on the beginning and ending dates; this practice was supported 
by inspection of contours of these values on maps of potential vorticity. 

Then, for each of the starting 0 surfaces, the parcels' ending 0 values were subtracted to 
get the average descent over this period. Figure 4 shows the result. Parcels starting at 600 K 


Figure 4 
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in mid-January had descended by some 40 K. on average, by mid-March. Parcels in the lower 
stratosphere, between 400 and 500 K, descended roughly 15 to 20 K during the same time. 
To get the average diabatic descent as a function of initial potential temperature, two straight 
lines were fit to the data, meeting at 475 K. 

To determine ozone loss, ozone fields were reconstructed on January 10 and March 17 
on a series of pairs of potential temperature surfaces within the vortex, the March 17 surface 
being adjusted upward to compensate for the descent of air within the vortex. The ozone 
decline was calculated by differencing these two vertical profiles; this method will be referred 
to as the “vortex average descent" method below. 

The second method used to remove diabatic effects to determine ozone loss involved 
tracing the parcel trajectories of the ozone measurements themselves. It the trajectory model 
were perfect, then parcels started on January 10 with certain PV and 0 coordinates would 
end on March 18 with different values of PV and 0. It these “diabatically adjusted PV-0 
coordinates are used in the linear time fits, then the slopes of those fits should leflect chemical 
loss alone. This method will be termed the “diabatic coordinate adjustment" technique in the 
discussion which follows. 

Of course, trajectories ot individual parcels cannot be traced accurately tor such long 
periods, but TRAJREF indicates that reasonable results can be obtained using large ensembles 
of trajectories whose results are averaged together. 
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4. Results 

4.1. Vortex Average Descent Method 

To apply the vortex average descent method of determining loss, ozone measurements 
from the ER-2 in situ instrument, the DC-8 AROTEL and DIAL lidar instruments, and the 
THESEO ozonesondes were combined into a single P V-0 composite. Random points from 
each instrument were discarded until their measurements were of approximately equal density 
in PV-0 space, so that no one instrument would dominate the analysis. The resulting composite 
was used to reconstruct ozone on potential temperature surfaces for January 10 and March 
17. The surfaces for March 17 were adjusted upwards to account for diabatic descent as 
determined by the fit in Figure 4. Reconstructed gridpoints of ozone in the core of the vortex 
were selected and area-averaged to obtain two vertical profiles, one at each end ot the time 
period being examined. 

Figure 5 shows these profiles. The difference between these curves, divided by the time 
period, is the average chemical ozone loss rate. The ozone loss rate calculated here peaks 
at 470 K with a value of 0.02G ppmv/day and decreased to near 0 at 650 to 700 K. Results 
are shown using the NCAR/NCEP Reanalysis; the UKMO analyses and DAO data produce 
similar curves with loss rates of 0.025 and 0.02G, respectively, at 470 K. 

The profiles in Figure 5 extend down to 350 K and show the loss rates continuing to 
shrink below the 400 K level. Although the PV-fl technique is more suitable for use in the 
stratosphere, where (modified) PV and 0 are at least quasi-conserved and quasi-orthogonal, the 
technique seeems to yield reasonable values down into the very lowest part of the stratosphere 


Figure 5 
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and into the uppermost troposphere. 

4.2. Diabatic Coordiate Adjustment Method 

To apply the diabatic coordinate adjustment technique, measurements from the THESEO 
sondes were used, since they were sparse enough that their trajectories could be traced easily. 
All days’ sonde profiles were run forward to March 18 to obtain their diabatically adjusted 
PV and 9 values on that date. They were also run backward to obtain their diabatically 
adjusted PV and 9 values on January 10. Parcels whose PV changed by more than 25% over 
this period were discarded to reduce the effects of parcels entering or leaving the vortex. 

The remaining measurements, with their adjusted P V-0 coordinates, were then used to 
construct a time-varying composite, and the slopes from the resulting time fits were obtained. 
Uncertainties in the slopes were also computed from the covariance matrices of the time fit 
coefficients (as in Meyer [1975]). 

The loss rates were mapped onto vertical profiles in the vortex core; Figure 6 shows the 
area-weighted average of these profiles. Figure 6a shows the results for the forward trajectory 
run on January 10. The loss rate peaks at 0.032 ± 0.007 ppmv/day at 460 K. This is larger 
than that computed from the vortex averaged descent method, but it peaks in about the same 
altitude. (Vortex averaged descent loss rates computed from the ozonesonde data alone had the 
same peak location and value as from the composite of all the instruments.) Figure 6b shows 
the results for the back trajectory run on March 17; the loss rate peaks around 0.30 ± 0.008 
ppmv/day at 440 K. The peak is about the same magnitude (indicating agreement between the 
forward and back trajectory runs), but it is located about 20 K lower in potential temperature 


Figure 6 
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because of diabatic descent. 

The back, trajectory loss rates mapped to the January 10 profile are very similar to the 
forward trajectory rates for that date. Likewise, the forward trajectory results mapped to March 
17 are similar to the back trajectory results for that date. That is, both trajectory runs show 
similar peak loss rates which descend in 0 over the time period examined. Figure 7 shows both Figure 7 
rates mapped onto February 13 profiles (the middle of the period). One difference between 
the two runs is that the forward trajectory profiles exhibit a small offset near the bottom of the 
profile (and are very nearly zero at the top), while the backward trajectory profiles show a 
small offset near the top of the profiles (and are very nearly zero at the bottom). 

Figure 8 shows the time fit slopes (from the backwards trajectory run) for the 440 K Figure 8 
surface as a function of equivalent latitude ( Butchart and Remsberg , 1986) for February 13. 

The area-weighted average loss rate within the inner vortex region on the 440 K surface is 
0.030 ± 0.009, as shown in the profile. Expanding the average to include most ot the vortex 

decreases the loss rate slightly to 0.028 ± 0.00 1 at 440 K. 

Results from the backwards trajectory run are similar within the vortex, but the change 
rates outside the vortex vary greatly between the two runs (from approximately -0.01 
ppmv/dat to +0.02 ppmv/day at 470 K). The near-zero loss rate shown in Figure 8 outside the 

vortex must be considered largely fortuitous. 

From Equation 1 , the loss rate of total ozone Or becomes: 

DOr/Ot = #0 vj \(!p)/0t 

= IF r b( PY. 0)<lp + \V x (Ps)0pJ0t 

Jo 


(4) 
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where U' = 0.789 m s 2 /kg and p s is the surface pressure. We do not reconstruct ozone values 
in the lower troposphere, but fixing the (notional) ozone profile to zero at the suiface (as is 
done to derive the total ozone values mapped in Figure 3) makes the second term vanish. (As 
a dynamical effect, the second term has no bearing on chemical ozone loss anyway.) Thus, 
integrating the time slopes yields that part of the time rate of change of total ozone which is 
caused by chemical destruction. 

Time fit slopes from the diabatically-adjusted sonde data were mapped onto vertical 
profiles as a function of pressure. Loss rates from the back trajectory results were used because 
they were close to zero below 400 K, a region more heavily weighted in the integration. Any 
ozone increases (i.e., from the offsets near the top of the profiles) were zeroed out— only ozone 
destruction rates are of interest. These rate profiles were then integrated and area-averaged 
within the vortex core for each day between January 10 and March 17. 

Figure 9 shows the daily total ozone loss rates for the core of the vortex. The losses 
begin at approximately 0.8 DU/day and increase to 1.3 DU/day by March 17. The loss rates’ 
evolution is associated with both a descent in the peak altitude of ozone loss and a shift in the 
PV values associated with the vortex core, with the latter dominating. Summing the daily loss 
rates yields a decrease of about 77 DU caused by chemical loss ot ozone. 

5. Discussion and Conclusions 

Ozone measurements from AROTEL, DIAL, the NOAA ER-2 instrument, and the 
THESEO ozonesondes have been combined using a PV-0 analysis technique which allows 


Figure 9 


the ozone field to evolve in time. Realistic ozone fields have been obtained by mapping the 
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ozone field from PV-0 coordinates back into real space. Diabatic effects have been estimated 
using modeled trajectories of large numbers of parcels. Estimates of chemical ozone loss in 
the lower stratosphere are found peak at 0.030 ± 0.009 ppmv/day near 460 K in January, with 
this peak descending to 440 K in March. The loss near 450 K is thus about 2 ppmv over the 
period. The loss rates also decrease to near 0 at 550 to 600 K. 

The P V-0 analysis has several potential pitfalls. It depends on its constituent being 
long-lived enough to be well-mixed along contours of isentropic potential vorticity. In the 
middle stratosphere and above, where photochemistry dominates, the method is of limited 
applicability to ozone. This analysis has been confined to the lower stratosphere. 

When an ozone measurement is made near the vortex edge in the sunlight, where 
chlorine-driven photochemical destruction is likely to be enhanced, this measurement may not 
be typical of ozone values around a PV contour on a 0 surface that intersects the measurement 
location. With enough additional measurements around that contour, PV-fl analysis will yield 
an averaged view of the ozone field which the ozone-depleted measurements will influence 
but not necessarily dominate. Poor sampling, however, may lead to the sunlit region being 
under-represented or over-represented in the average until the ozone-depleted air is mixed 
along the contour. This analysis relies on the time series fit to average out the sudden drops in 
ozone that might appear in such cases. 

To distinguish between diabatic effects and chemical destruction of ozone, this analysis 
uses diabatic trajectories of large numbers of parcels traced over about ten weeks. Averaging 
the effects of large numbers of parcels might well compensate for the inaccuracy of 
individual parcel trajectories, but there remains the possibility of systematic artifacts ot the 



18 


meteorological analysis used or of the radiation model used to generate the heating rates. 

The diabatic heating rates have been verified in previous work ( Rosenfield et al. 
[submitted]) for the UKMO meteorological products, and they appear to be reasonable. 

Another approach would be to analyze a tracer such as N 2 0 and test whether the diabatic 
trajectories reduce the apparent rate of change, which should be a purely diabatic effect. This 
has been attempted, and the results suggest that using diabatically-adjusted trajectories tends 
to reduce the rates of change in the vortex core at altitudes higher than 500 K. However, the 
N 2 0 measurements from SOLVE were sparse enough, and the uncertainties in the rates are 
consequently large enough to make these results inconclusive. That is, the differences between 
the uncorrected rate profiles and the diabatically corrected rate profiles are smaller than the 
uncertainties in either. 

Another effect which would complicate matters is mixing across the vortex boundary. 
Although the polar vortex is generally recognized as being isolated from the midlatitudes 
(Bowman [1993], Schoeberl et al. [submitted]), it is not completely isolated— slow or weak 
mixing may occur over the course of a season, in addition to occasional major intrusion events 
such as one described by Plumb et al. [1994], The PV-tf technique may still reconstruct ozone 
fields with some degree of fidelity in these circumstances, but untangling chemical destruction 
from mixing becomes problematic. In this analysis, the issue is addressed by excluding 
parcels whose diabatic trajectories involved a change of PV of more than 25%, but it is unclear 
whether this is sufficient given the limitations of tracing individual parcels. In other words, 
it is difficult to say whether a parcel has been excluded because it moved into or out of the 
vortex, or because its trajectory is simply not accuiately traced. 
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Inspection of the gridpoint trajectories used to determine descent rates shows only a 
few cases where trajectories which originated outside the vortex in January ended up inside 
the vortex in March: about 1.3% of the mid-March vortex parcels at 420-430 K for the 
back-trajectories only. For other theta levels and for the forward trajectory run, no such 
parcels were found. Schoeberl etal. [in press], which used a more comprehensive trajectory 
modelling run, supports the assertion that mixing into the vortex was negligible over most of 
this period. 

Toward the middle of March, as the vortex was breaking up, the effects of mixing are 
expected to increase. A sharper drop-off in ozone at the end of the time series can be seen in 
Figure 1, and it is difficult to say how much of this is caused by mixing and how much by 
chemical destruction in the increased sunlight. This analysis partially side-steps the issue by 
fitting a linear trend over the 10- week period, so that the sudden drop-off at the end does not 
influence the overall trend too strongly. 

Other techniques have been used to categorize ozone loss. Instead of placing all the 
data into a composite time-dependent field, the MATCH program (Rex et al . , 1999) uses 
paired individual ozonesonde profiles to determine ozone loss. MATCH selects the pairs by 
tracing the trajectories of sonde-sampled air parcels until they intersect other sonde launches 
days later. By using small clusters of trajectories for each measurement, and by averaging 
differenced pairs together, MATCH tends to reduce the uncertainties associated with the 
trajectory calculations. Because MATCH uses changes in PV as a criterion in matching pairs 
of sonde launches, its ozone values are effectively categorized by PV in a way similar to PV-0 
analysis. Instead of dealing with pairwise differences at discrete 0 levels, though, the P V-0 
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technique has the effect of averaging all ozone measurments of similar PV and 0 values. More 
data will tend to enter such an average at the cost of a higher dependence upon the assumption 
than ozone is well-mixed along an isentropic PV contour. 

The vortex-averaged analysis of Schoeberl et al. [submitted] involves diabatic trajectory 
calculations initialized from massive numbers of ozone measurement locations, so that the 
inaccuracies associated with individual parcels are averaged out. The data which remain 
inside the vortex are fit to a quadratic curve at various theta levels to obtain ozone loss rates 
from December 1999 through mid-March 2000. Solar exposure is also used to discriminate 
between parcels in the core of the vortex and parcels from the edge region. 

The PV-0 technique with diabatic corrections may be thought of as a way ol enhancing 
such a trajectory analysis to produce a continuously evolving ozone field. The peak ozone loss 
rates produced here are lower than MATCH’S of 0.05 ppmv/day, but somewhat higher than the 
0.025 ppmv/day obtained by the Schoeberl trajectory analysis. The peak loss rates seem to be 
a little lower and broader in altitude in the latter results, as well, but this may be an artifact ot 
the heavy smoothing done here. 
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Figure 1. Time series of ozonesonde mixing ratio measurements near PV=32 PVU and 
0=450 K. The gray shading of each symbol indicates its weighting in the linear time fit; the 
size of each symbol indicates its proximity to the P V-0 gridpoint. The linear time fit and its 
uncertainty (± one standard deviation) are overlaid. Dashed vertical lines mark February 1 and 
March 1 . 

Figure 2. Ozone as a function of flight time and altitude for the DC-8 flight of 2001-03-09, (a) 
as measured by the AROTEL instrument on board the DC-8, (b) using THESEO ozonesonde 
measurements mapped into PV-0 space and then mapped onto the AROTEL measurement lo- 
cations. 

Figure 3. (a) Total ozone fields as measured by TOMS (left) and as calculated from sonde, 
lidar, and ER-2 ozone data using the P V-8 reocnstruction technique (right), (b) Scatter plot of 
reconstructed values versus the TOMS measurements. 

Figure 4. Vortex descent between 2000-01-10 and 2000-03-18 as a function of starting poten- 
tial temperature. Crosses mark parcels from forward trajectories, and diamonds mark parcels 

from backward trajectories. Two linear least-squares fits have been drawn through the points. 
Figure 5. (a) Reconstructed ozone profiles (area-averaged over the core of the vortex) from 

January 12 (right line) and March 17 (left line). The March profile has had its 0 values adjusted 

to remove diabatic descent over the period, (b) Ozone loss rates for the voitex core, calculated 

by subtracting the two curves from (a). 

Figure 6. Ozone loss rate (from the time fit slopes of diabatically-adjusted sonde data) as 
a function of 0 for the vortex core region for (a) January 10 using diabatic corrections from 
forward trajectories and (b) March 17 using diabatic corrections from back trajectories. 
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Figure 7. Same as Figure 6, except both ozone loss rates have been mapped for February 13 
for both forward and back trajectory results. 

Figure 8. Ozone loss rate (from the time fit slopes of diabatically-adjusted sonde data) as a 
function of equilalent latitude at the 440 K 9 surface on February 13, 2000. The vertical dotted 
line marks the core of the polar vortex. 

Figure 9. Vertically integrated chemical ozone loss rates for January 10 through March 17, 
within the core of the polar vortex. 
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